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RECURRENCE RELATIONS FOR EXCEPTIONAL HERMITE 

POLYNOMIALS 


D. GOMEZ-ULLATE, A. KASMAN, A.B.J. KUIJLAARS, AND R. MILSON 


Abstract. The bispectral anti-isomorphism is applied to differential opera¬ 
tors involving elements of the stabilizer ring to produce explicit formulas for 
all difference operators having any of the Hermite exceptional orthogonal poly¬ 
nomials as eigenfunctions with eigenvalues that are polynomials in x. 


1. Introduction 

The classical orthogonal polynomial families of Hermite, Laguerre and Jacobi m 
have three important properties: they are eigenfunctions of a differential operator 
in X having eigenvalues dependent on n, they satisfy a well-known three-term re¬ 
currence relation, and they form a basis of a given space. 

In the past few years, the class of polynomials having the above properties has 
been extended to include exceptional orthogonal polynomials, which differ from 
the classical families in that there are a finite number of degrees for which no poly¬ 
nomial eigenfunction exists. The differential equation for exceptional polynomials 
contains rational instead of polynomial coefficients. If certain regularity assump¬ 
tions are imposed then, like their classical counterparts, the exceptional orthogonal 
polynomials form the basis of a weighted space. 

This paper is mostly concerned with the study of recurrence relations for ex¬ 
ceptional Hermite polynomials. Some particular examples of these polynomials 
were investigated as early as [3], but a final classification appeared only recently 
in [3] . See that paper for a more detailed overview of the subject and for additional 
references. Exceptional Hermite polynomials were also derived from exceptional 
Charlier polynomials by taking a suitable limit [5]. 

Exceptional polynomials admit two different sets of recurrence relations. The 
first set of recurrence relations MM have coefficients that are polynomial functions 
of X and n and consequently have no obvious bispectral interpretation; i.e., they 
cannot be interpreted as higher-order Jacobi operators. We shall not investigate 
these relations in this paper. The second set of recurrence relations have order 
2N + 3 where N is the codimension (number of missing degrees) of the exceptional 
family. Except for one term, the coefficients do not depend on x and are rational 
functions of n. This type of relations will be the main focus of our study. 

The existence of this second type of recurrence relations of order 47V -|- 1 was 
first established in ITO]. Relations of order 27V -|- 3 order for exceptional families 
obtained through a one-step Darboux transformation were recently given in m, 
and later generalized in but in both cases an explicit construction of the 
recurrence coefficients was missing. Somewhat closer to this goal is the recent work 
of Duran [7], which provides explicit relations for some examples of exceptional 
Hermite and Laguerre polynomials. His approach uses a bispectral correspondence 
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between discrete Krall polynomials and exceptional Charlier polynomials and a 
limit procedure to obtain the recurrence relations for the exceptional Hermite. 

The main result of this note (Theorem 15.81) is a straight-forward procedure for 
producing the general recurrence relations satisfied by the exceptional polynomials. 
We focus on the Hermite subclass but it is clear that the same techniques can be 
applied to the other families as well. In particular, by making use of a bispectral 
anti-isomorphism of rings of operators acting on the classical Hermite polynomials, 
the production of recurrence relations for the exceptional Hermite polynomials is 
reduced to elementary algebra. 

The idea of the bispectral anti-isomorphism for creating “bispectral Darboux 
transformations” was pioneered by Kasman-Rothstein [13] . further developed by 
Bakalov-Horozov-Yakimov [2], and first used in the context of orthogonal polyno¬ 
mials by Griinbaum-Yakimov The construction below does not depend on 

those prior results but is an application of the same ideas to the new context of 
exceptional orthogonal polynomials. 


2. Background 


2.1. Hermite Polynomials. Set h{n,x) = Hn{x), where the latter is the nth 
degree Hermite polynomial as defined in Section 18.3 of [18]. These polynomials 
have associated operators with the following properties. 

Letting 

= f{n + l) 

denote the elementary shift operator, define the following degree-raising and lower 
operators, 

(2.1) A = i0-hno0-i, 

(2.2) r = 2no0-i, 

where the n denotes the indicated multiplication operator. The action of A and 
r on the classical Hermite polynomials can also be expressed as operators in x, 
namely 

(2.3) xh{n,x) = Ah{n,x) 

(2.4) dh(n,x) = rh{n,x) 

where d = -^ differentiates with respect to x. Note that (12.31) is just the 3-term 
recurrence relation for the classical Hermite polynomials. We will also need the 
second-order operator 

T{y{x)) = y"{x) - 2xy\x) 
and express the Hermite differential equation as 

(2.5) T{h{n,x)) =—2nh{n,x). 

Remark 2.1. To read and understand the remainder of the paper, it is helpful to 
keep in mind that operators in x only always commute with operators acting only 
in n. In general, we will try to write operators in x with Roman symbols (except for 
d) and operators in n with Greek ones to help the reader remember which variable 
an operator acts in. However, the map b to be introduced later turns operators in 
X into operators in n. Moreover, b(i?) and R have the same effect on h{n,x) for 
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any R G <C\x,d]. So, when either one of these expressions appears next to h(ri,x) 
in a formula it is possible to replace it by the other. 

2.2. Exceptional Hermite Polynomials. This subsection quickly reviews the 
definitions related to the Exceptional Hermite Polynomials h{n, x) associated to 
the choice of a partition A. For proofs of the claims and more details, see [9]. 

Let Ai < A 2 < ■ • • < A^ be a non-decreasing sequence of positive integers so that 
A = {Ai,..., Af} is a partition of the positive integer N = into i parti3- It is 
convenient to associate to the choice of A the set K = {ki,, ki}, where 

(2.6) ki = Xi+i — l, i 

Using Wr to denote the Wronskian determinant with respect to x, to the selection 
of partition A (or sequence K) we associate the polynomial r]{x) defined by 

(2.7) r]{x) = Wr [h{ki,x),..., h{ki, x )], 

and define the order £ differential operator A by its action on an arbitrary function 

y{x), 

(2.8) = Wr [h{ki,x),..., h{kt, x), y{x)] . 

Now, a new sequence of polynomials is defined by 

(2.9) h{n, x) := A(h(n + £ — N,x)), 

where it is clear from ()2.8I) that h(n, x) is either a degree n polynomial in x or zero. 
The degree sequence of the resulting non-zero polynomials is 

(2.10) I = {nGNo:n>N-£,n + £-N^K}. 

The degrees 0, l,...,iV —£—1 and the degrees ki+N — £,...,ki + N —£ are missing 
from X, so that the polynomial sequence {/i(n,x) : n G X} is missing a total of N 
degrees. Let us also denote hy 14 C C[x] the set of finite linear combinations of 
exceptional Hermite polynomials: 

(2.11) U = spaii{h(n,x): n G X} 

The resulting families of polynomials generalize the classical Hermite family be¬ 
cause they satisfy a Hermite-like differential equation 

(2.12) T{h(n,x)) = 2{N — n)h(n,x), nGX 
with 

(2.13) fiy) := y" - 2 fx + j/' + f ^ + 2x^') y. 

\ V J \V V J 

We will say that A is an even partition if £ is even and if A 2 i-i = X 2 i for every 
i < £/2. For an even partition A, the exceptional Hermite polynomials h(n,x) 
satisfy the orthogonality relationtO 

r ^ ^ ^ ^ 

/ h(n,x)h{m,x)W{x)dx = Sm ~ ki), n,mGX 

•IR ’ i=i 


^Many of the objects in this paper depend on the selection of A. However, to avoid cumbersome 
notation, we will consider A to be selected and fixed and not specifically indicate the dependence 
of objects, such as h{n, x), S, U and X on this choice. 

^This formula corrects a misprint in equation (52) of [^. 
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where j = n + £ — N, and the weight is defined as W{x) = j■ Moreover, 
if A is an even partition then TA is dense in W{x)dx). 

Despite the importance of even partitions to the role of h{n, x) as orthogonal 
polynomials, the formulas for the recurrence relations that constitute the main 
result of this paper are valid for an arbitrary choice of A so this distinction will not 
be further emphasized. 

2.3. Bispectrality. The bispectral problem, as proposed by Griinbaum and ini¬ 
tially studied by Duistermaat-Griinbaum [5] seeks to find functions depending on 
two variables that are joint eigenfunctions for operators in each of the two variables 
separately so that in each case the eigenvalue depends non-trivially on the other 
variable. Specifically, a bispectral triple {L,A,ip) consists of a function ip(x,z), 
an operator L in x satisfying L'0 = p{z)'tp and an operator A in z satisfying 
Kip = TT{x)ljj. 

A great body of research by many authors has investigated bispectrality for vari¬ 
ous different kinds of operators, establishing theorems, classifying examples, finding 
connections to other areas of research such as soliton equations, or quantum and 
classical particle systems. (See, for instance, [H] and references therein.) Addi¬ 
tionally, it has long been recognized that the classical orthogonal polynomials can 
be interpreted as exhibiting bispectrality if they are viewed as functions of two 
variables (one continuous and one discrete) rather than as a sequence of separate 
functions of x. The differential equation in x and the recurrence relation are then 
the two eigenvalue equations. (See, for example, [TUI 1^ .1 

3. The Anti-Isomorphism 

The key observation leading to the main result is that the ring C[a;,9] of poly¬ 
nomial coefficient differential operators and the ring C[A,r] of difference operators 
generated by A and T are naturally anti-isomorphic. 

Definition 3.1. Define the map b by its action on monomials: 

\>:C[x,d] C[A,r] 
x^d^ hA PA* 

extended linearly to all of C[x, 9]. 

Theorem 3.2. The map b : C[x,9] —>■ C[A,r] satisfies 

Q{h{n,x)) = \>{Q)h{n,x) VQ S <C[x,d\ 

and 

b(Q o i?) = b(i?) o b(Q) VQ, R € C[a:, 9]. 

In other words, b is an anti-isomorphism that has the property that an operator and 
its image always have exactly the same effect on the classical Hermite polynomials. 

Proof. Because b extends linearly it is sufficient to check the first property on a 
monomial of the form x'^dP Then observe that because x and T commute one has 

x^d^h{n,x) = x*P7i(n,x) = x^h{n,x) = lK^h{n,x) = \>{x^ d^)h{n,x). 

As a consequence, we have for any Q and R in C[a:, 9] that 

b(Q o R)h{n, x) = {Q o R)h{n, x). 
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But, note that 

b{R) o b{Q)h{n, x) = \){R){Qh) = Q{}){R)h{n, x)) = Q{R{h{n, x))) = {Q o R)h{n, x) 

also because Q commutes with b(i?). Then b{Q o R) — b{R) o \){Q) is an operator in 
C[A, r] which has the function h(n, x) in its kernel for all n. This implies that it is 
the zero operator. □ 


4. Factorization Lemma and Eigenvalue Equations 


Define another differential operator B of order ^ by stating that for an arbitrary 
function 2 /(x), 


(4.1) 


B{y{x)) = (-1) 


r;“^Wr 


hi{x),... ,he{x),e ^ y{x) 


where 

(4.2) hj { x ) :=Wi 

with the hat decoration indicating omission from the sequenc^. The main result of 
this section is Corollary 14.41 providing eigenvalue equations for the products Ao B 
and B o A. Unfortunately, the introduction of several technical definitions and 
lemmas is necessary to obtain it. 

For 0 < j < k, let denote the truncated partition = (Ai < • • • < Aj) 
and = {fci,..., fcj} the corresponding truncated set. Let 


h{ki,x),... h{kj,x), ..., h{kg, x) 


(4.3) 


r]j{x) =Wi[h{ki,x),... ,h{kj,x)] , 


denote the corresponding Hermite Wronskians. For 1 < j < i, define the first-order 
differential operators 


(4.4) 


(4.5) 


Aj{y) ■= 
Bj{y) ■■= 


Vj 

Vj-i 

Vj-i 

Vj 




Wr [r]j,y] 


Vj-i 


y' - (2x + 


Vj-i 


Wr 


y = 


e Vj-i,y 


e-^ rjj 


where we define rjo = 


Proposition 4.1. The above operators are linear factors of the order ^ operators 
A and B defined in (US) and dsi). Indeed, the following factorizations hold: 

(4.6) A = Ai o ■ ■ ■ o A 2 o Ai 

(4.7) B = Bi o B 2 o ■ ■ ■ o Bi 


Proof. The proof makes use of the following classical identities m- Throughout, 
the arguments are sufficiently differentiable functions of one variable. 


(4.8) 

(4.9) 


Wr [g/i, 5 / 2 ,..., gfn] = 5 " Wr [/i, / 2 , 


Wr [Wr [/i, 5 ], Wr [/i,..., /„, h]] = Wr [/i,..., /„] Wr [A, ..., /„, g, h ], 


(4.10) 


Wr 


/i,...,/J =(-l)L-/2JWr [A,..., 


^This formula corrects a misprint in equation (27) of [^. 
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where 


fj = Wr 


Let us prove (14.61) . Since Aj annihilates rij and Aj-i o ■ ■ ■ o Aih(kj,x) = r]j{x) 
(which can be shown using (14.3L (14.41) and (I4.9p l. the product • -oAi annihilates 

all h{kj,x), j = 1,... ,£. By construction, the £th order operator A has the same 
property and so must be a multiple of the product. Equality is established by 
considering the coefficient of the highest order derivative y^^\x). 

We prove (IQ) by induction on For £ = 1 the equality of (HU) and (1431) 
with j = 1 follows from (14.8|) . Suppose that the identity (14.71) holds for £ — 1. For 
1 < * < 4 < f'l let hij denote the £ — 2 order Wronskian like (14.21) . but with h{ki, x) 
and h(kj,x) omitted. Since 

Tlt-i = hi, 

we have by (14.51) . (14.81) and (14.91) . 


Bi{e^ hj) = 


■ Wr 


m 


hi, hj 


= -e^ hi 


By the inductive hypothesis, 
{Bio---o Bi^i) 


hi,i, • • • , hi-i^i, hjj 


= 0 , 


j = l,...,£-l, 


and hence 


(Bi o • • • o Bi^i o Bi) 


= 0, j = 1, 


The last identity also holds for j = £, since Bi{e^ hi) = 0. 
By (|4.1I) we also have 


B 


(e^^h,) 


= 0, j = l,...,J 


Thus Bi o ■ ■ ■ Bi and B are two £-th order linear differential operators that agree 
for £ linearly independent functions and hence they agree up to a multiplicative left 
factor. Directly from the definition (14.51) . we have 

[Bi o ■ ■ ■ o Bi){y) = — + terms with lower order derivatives. 

m 

By (14.11) . (14.81) and (14.101) . we have 




Wr 


B{y) = (-1) 


hi,... ,e^ hi,y 


= (_1)L^/2J. 

1 


Wr 


(e^'^rji) 
hi,..., hi 


M) 


iVi) 


terms with lower order derivatives 


= —y 
m 


w 


terms with lower order derivatives. 


This establishes (1431) . 


□ 
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Next, let 

T,(y) := y" - 2 (x + y' + (^ + y, 

\ vjJ \ Vj m J 

denote the exceptional operators corresponding to the truncated partition 
Thus To = T is the classical Hermite operator, while Ti = T, see (I2.13L 

Proposition 4.2. With Aj,Bj,Tj defined as above, the following intertwining re¬ 
lations hold: 


(4.11) Ai,oT,_i = (T,-2)oAi„ 

(4.12) T,_i o B, = B, o (T, - 2). 

Proof. Forj = 1, ...,£, set (a;) = x) and observe that Wr [i/^i(a;),...,'^^(a:)] 

Next, set Hj = —d^ + Uj, where 

Uj = x^ - 2dl log rij{x) + 2j. 

The operator Bj is such that 


(4.13) 

Since 


„-x^/2 


Vj 


O Bn O 


,-x^/2 

— 


— —Tj + 2_7 + 3. 


{-dl + x^)f)j = {2kj + l)f}j 


Crum’s Theorem [3] implies that 

Wr[V>i,.. 


B,-i 


Equivalently, in view of (14.311 . 


— {2k j + 1) 


Wr[V>i,...,i/;j-] 


B,-i 

This gives the factorizations 
Bj-i = i-dx + a; - 


rij-i 


— ( 2 % + 1 ) 


B-i 


Vj Vj-i 


o \dx + x - - + ■ 


Vj V]-i 


Bn = [dx-kx - - 


% , ^a-i 


Vj Vj-i 


° -dx + x - - 


^j , ^j-i 


Vj Vj-i 


+ 2kn 4“ 1, 


4“ 2kn 4“ 1. 


Hence Bj is obtained from Bj-i by a Darboux transformation. Straightforward 
calculations show that 


^j — 


Vj 


Vj-i 
I Vj 


dx-^ 
Vj 
-1 


oldx+X-"^ I 


Vj Vj-i 


-Bj = 


!M^(-dx + 2x+'^) 
Vj V Vj-iJ 


e^ 

Vj-i 


-1 


o -dx + a;- - 


^j , ^j-i 


Vj Vj-i 


Vj-i 


— 
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Thus, for 1 < j < i, 

(4.14) T,_^=B,oA,+2{j-k,-l), 

(4.15) Tj = o + 2(j - kj). 

Multiplying (j4.14|) on the left and (14.151) on the right by Aj and equating the terms 
equal to Aj o Bj o Aj yields the desired equalities. □ 

The operators A o B and B o A are polynomials in the operators T and T 
respectively: 

Lemma 4.3. The products of the operators A and B defined in (|2.8I) and 
factor as: 

(4.16) BoA={T + 2ki)o---o{T + 2ke), 

(4.17) AoB = {f-2e+ 2ki) o---o{f-2i + 2ke). 

Proof. It follows from (|4.14l) that 

B\ o Ai = Tq — 2ki = T — 2ki. 

Now, suppose that for some 1 < n < ^ that 

(4.18) BioB 20 - ■ -oBnOAnOAn-lO- ■ ■ O Ai = (T +2ki)o (T + 2 ^ 2 ) 0 - ■ -oiT +2kn). 

Consider the left-hand side of (14.181) with n replaced by n-l-1. Using (14.141) once to 
replace Bn+ioAn+i and following repeated use of (14.111) to move the term involving 
Ti all the way to the right we obtain 

Bi O ■ ■ ■ O BnO {Bn+l O A„+i) O O • • • O Ai 
= Bi O ■ ■ ■ O BnO {Tn - 2(n - kn+l)) o A„ O • • • O Ai 
= Bi O ■ ■ ■ O BnO An o {Tn-1 - 2{n - kn+l - 1)) O • • • O Ai 


= Bi O ■■■ O Bn O An O ■■■ O Ai O {Tq + 2kn+l) 

= (T -\- 2k\) o {T -\- 2 ^ 2 ) o • • • o {T -\- 2fcn-i-i) 

where the last equality is obtained by making use of the assumption (14.181) and the 
fact that Tq = T. 

Then by induction we find that (14.181) is true when n = £. Since the left-hand 
side in that case is equal to o A by Proposition 14.11 this proves (14.161) . 

Factorization (|4.17|) can be similarly obtained by repeated use of (|4.12|) and 

( 144 ^ . □ 

Combining Lemma [4.31 with the eigenvalue equations p.5|l and (I2.12|) leads im¬ 
mediately to: 

Corollary 4.4. The operators Ao B and B o A satisfy the eigenvalue equations: 

(4.19) {B o A)h{n,x) = TT{n)h{n,x), 

(4.20) (A o B)h{n, x) = 7r(n — iV -|- £)h{n, x). 

where 

t 

’^{n) = ]^(~2n -I- 2ki). 

i=l 


(4.21) 
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5. The Stabilizer and Recurrence Relations 
Definition 5.1. Let us define the stabilizer ring 5 C C[a;] by 
5 = |/ G <C[x\ : f{x)h{n,x) G U for all . 

Remark 5.2. By the orthogonality properties of the polynomials h{n, x) it is possible 
to show that if / G 5, then for all n G I, necessarily 

deg/ 

f{x)h(n,x)= ^ jj^riHn + j,x) 

j=- deg/ 

for some G C. 

It was already observed in m that rj^{x) G S and more recently in [Biiaiii] 
that 77 G 5 as well. The following proposition provides a characterization of S. 

Proposition 5.3. If f{x) G C[a;] and f'{x) is divisible by f]{x), then f G S. The 
converse is also true, provided rj{x) has only simple roots. 

Proof. We first show that p G C[a;] belongs to U if and only if 

(5.1) 2r]'{x){xp{x) —p'{x)) + p" {x)p{x) 

is divisible by p. Indeed, the operator (I2.13P transforms every element p G U into a 
polynomial, so the numerator of the singular part of T(p) must be divisible by p: this 
is precisely (EB. Conversely, divisibility of m by p imposes precisely deg p = N 
independent linear conditions on p G C[a;]. Since N is also the codimension of U it 
follows by dimensional exhaustion that a polynomial p G C[a;] that satisfies these 
conditions is necessarily an element of U. 

Now suppose that f'{x) is divisible by p{x). Let p GU and set q = fp. Then, 

2p'{x){xq{x) - q'{x)) + p"{x)q{x) = 

f{x){2p'{x){xp{x) - p'{x)) + p"{x)p{x)) - 2p'{x)p{x)f'{x) 
is divisible by p, and hence is also in U. 

Conversely, suppose that p{x) has simple zeros, and that the above expression 
is divisible by p{x) for all p gU. Also, p{x)f'{x) is divisible by p{x) for all p GU. 
By Proposition 5.5 in [9], the polynomials in Li do not have a common root, which 
implies that f'[x) must be divisible by p{x). □ 

Note that as a consequence of Proposition 15.31 the ring S is not only non-empty 
but moreover guaranteed to have elements of every sufficiently high degree. Hence 
if p{x) has simple zeros, then f^pGSis the element in S of minimal degree]^ 
The elements of S can also be characterized as precisely those polynomials in 
C[a:] that act as linear endomorphisms on U. From we have 

Lemma 5.4. The operators A,B defined in (ESI) and m act as linear transfor¬ 
mations A : <C.[x\ ^lA and B : U —>■ C[a;], respectively. 

An immediate consequence is the following. 

^The minimal degree of an element of S and thus the minimal order of the recurrence relation for 
the exceptional Hermite polynomials was conjectured to be 2AI-I-3 in the previous works 171141171 . 
We see that this conjecture on the minimal order relies on Veselov’s conjecture on the simplicity 
of the zeros of the Wronskian of Hermite polynomials [8] . 
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Lemma 5.5. For fGS,BofoAisa differential operator in C[a;, d\. 

Proof. By definition of S and the preceding lemma, applying B o f o A to any 
polynomial results in a polynomial. Applying it to 1 then reveals that the order 
zero term is a polynomial. The claim then follows by induction since applying the 
operator to reveals that the coefficient of is a linear combination of powers of x 
multiplied by the lower order coefficients and the result, all of which are polynomial 
by assumption. □ 

The reason we need the previous lemma is so that we can apply the anti¬ 
isomorphism b to B o f o A: 

Definition 5.6. For f G S, let 

(5.2) A/ = 0^“^ob(Bo/oA)o-^o0^-^ 

Lemma 5.7. The linear transformation B : U ^ <C[x\ is injective. 

Proof. By (|4.20p . 

(A o B)h(ri, x) = 2^ (n(t. + N — t — n)^ h{n,x). 

By (12.101) . if n € I is one of the permitted degrees, then the right side is not zero. 
Hence, none of the basis elements of 14 can be in the kernel of B. □ 

Our main result is that the differential operator Af has h(n,x) as eigenfunction 
with eigenvalue /: 

Theorem 5.8. For all f G S, Afh{n,x) = f{x)h{n,x). 

Proof. Let f G S and set 

A/ = b(H o/o A) o 

7r(n) 

Hence, applying A/ o 7r(n) to h{n,x) and invoking Corollary 14.41 sfives 
{Af o Tr{n)){h{n, x)) = {Af o B o A){h{n, x)) 

= {B o Af o A){h{n, x)) 

(5.3) = {B o Af){h{n + N — i,x)) 

Furthermore, 

b(i? o / o A){h{n, x)) = {B o f o A){h{n, x)) 

(5.4) = {B o f){h{n + N — £,x)) 

Since A/ o 7r(n) = \>{B o f o A), the expressions (15.3p and (15.41) are equal and so the 
operator B annihilates the polynomial 

Af{h{n + N — x)) — f{x)h{n + N — £, x). 

For each fixed n, the first term is a linear combination of polynomials in lA and 
the second term is in W, by assumption. Since no nonzero elements of lA are in the 
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kernel of B (see Lemma [5.7L this must be the zero polynomial. Finally, note that 
Afh{n,x) = o Af o Q^~^)h{n,x) = {Afh{n + N — £,x)) 

_ (^f(^x)h{n + N — £, x)) = f{x)h{n, x), 

as claimed in Theorem 15.81 □ 

6 . Examples 

In this section we illustrate Theorem 15.81 by exhibiting a closed form description 
of the recurrence relation including some specific examples as well as some very 
general ones. 

6.1. One-step examples. Let us first consider the case of £ = 1 and write K = {k} 
where fc > 1 is an integer. By Proposition l5.3l and the fact that = 2{k + l)Hk, 
we have that h{k + 1 , cc) G 5, which leads us to the following identity. 

Proposition 6.1. Fix k > 1 and write 

h{n^ x) = Wr [h{k, x), h{n — fc + 1, cc)] , n>k — n^2k—1 

for the corresponding 1-step exceptional Hermite polynomial of degree n. Then, 

( 6 . 1 ) h{k l,x)h{n,x) 

^+1 /L I i\ 

= (n — 2/c + 1) 2-^ ( , j {n3 — k — j)j_ih{nk1 — 2j,x) 

1=0 k / 

where 

{X)n=X{X + l)---{X + n-l) 
denotes the usual Pochhammer symbol. 

The proof of Proposition 16.11 requires the following two lemmas. 

Lemma 6.2. Setting f{x) = hik + l,x) and 

= {y' - 2xy)/h{k,x) 

My{x)) = Wr [h{k, x), y{x)] 
as per the definitions above, we have 

[Bo f o A)[y\ = h(k -\-1, x)y" — h{k -\-2,x)y' -\-2k{h{k-\-\,x) — 2{k-\-l)h{k — l,x))y 
Proof. A direct calculation suffices to establish this. □ 

Lemma 6.3. For integers k,n>0 we have 

hik, A) = ^ 2^ n in-j + 1), o ©'=-2^ 

1=0 

where (n — j + l)j is the multiplication operator by the indicated polynomial in n. 

Proof. This is just the restatement of the linearization formula for the product of 
Hermite polynomials [1] p. 42] □ 


12 


D. GOMEZ-ULLATE, A. KASMAN, A.B.J. KUIJLAARS, AND R. MILSON 


Proof of Proposition \6.1\. By the first lemma, 

\){B o / o a 1) = o h{k + 1 , A) — r o h{k + 2, A) + 2k{h(k + 1, A) — 2(k + \)h(k — 1, 


Using the second lemma and elementary algebraic manipulations gives 

fc +1 /7 _|_ 1 \ 

\){B o f o A) = {k — n) ^ — 2j + 1) f . j{n — j + 

j=o \ J J 


Now 

7r(n) = 2(k — n). 

Applying the Definition 15.61 amounts to division of the preceding expression by 
—2{n — 2j + 1 ) and gives the desired formula. □ 


6.2. Multi-step examples. Next, we rederive the Miki-Tsujimoto-Duran relation 
[THU]. In this case, A = (1,1) and K = {1,2}, and 

(6.2) 7 r(n) = 4(n — l)(n — 2). 

We have 


(6.3) ri{x) = Wr [h(l, x), h{2, x)] = 4(1 -|- 2x^) 

(6.4) hiji, x) = Wr [/i(l, x), h{2, x), h(n, x)] . 

In order to avoid the use of fractions, we take 


o /-x 

(6.5) f{x) — - / r]{s)ds = x{3 + 2x'^) € S. 

4 Jo 

Applying (12.81) and (14.11) gives 

My) = 4((1 + 2 x^) 2 /" - 4 x 1 /' + 4y) 

4(1-h 2x2)^ (1-P 2x2)2 ^ 2 

Further calculations using this and (16.51) then show that 


( 6 . 6 ) {B o f o A){y) = x(3 + 2x^) y"” + (6 - Sx'^) y'" 

— x (6 + 8 x^ — 8 x'^) y" + x^(24 — 16x^) y' — x(24 — 16x^) y 

which has indeed only polynomial coefficients, as claimed in Lemma 15.51 
Then applying b according to our Definition 13.11 we find 


b(B o / o A) = r'‘A(3 + 2A2) + - SA'^) 

- r2A(6A + SA^ - 8A'‘) + 8rA^(24 - 16A^) - A(24 - IdA^). 

We rewrite this in terms of powers of the shift operator 0 by means of the identities 
(lO) (lO). This is a lengthy calculation (which is easier to do with the help of 
symbolic software) and it results in 

(6.7) 

b(B o / o A) = (n - 2)20^ + 6(n - 2)30 + 12(n - 3 ) 40 ”^ + 8(n - 5)2(n - 2 ) 30 - 3 . 


Note that we use Pochhammer symbols. 
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( 6 . 8 ) 


(6.9) 


By Definition 15.61 and (16.2p we have (since iV = ^ = 2 in this example) 

1 


4A f = [>{B o f o A) o 


{n - 2)2 


2)2 q3 _ 1 ) 20-1 + 8(n - 2 ) 30-3 


(n + 1)2 

where we used 0^ o a{n) = a(n + k) 0* and simplifications of the Pochhammer 
symbols. We therefore recover the desired recurrence relation 

(6.10) 4x(3 + 2x^)/i(n, x) = 4A//i(n, a;) 

(n - 2)2 : 


(n + 1)2 


h{n + 3, x) + 6(n — 2)h{n + 1, x) 

+ 12(n — l) 2 ^(n — l,a;) + 8(n — 2 ) 3 ( 1(71 — 3, x). 


The above procedure is entirely algorithmic and can be easily implemented using 
a computer algebra system. Here, for example, is the 2-step relation corresponding 
to A = (2,2). Let 

h{n,x) = Wr [/i(2, x),h{3, x), h{n — 2, x)] 
be the corresponding exceptional Hermite polynomial of degree n. Then 


32 


24a; + —x° I h{n, x) = - 


1 (n - 5)2 


5 ( 71)2 


h{n -I- 5, a:) -|- 2 


{n - 5)2 

71 — 1 


h(n + 3, x) 


+ 8(n — 5)(n — 3)h(n -b 1, a;) -b 16(71 — 2)^(77 — 4)/i(n — 1, x) 

32 

-b 16(71 — 5 ) 4 / 1(77 — 3,a;) -b —(77 — &)^h{n — 5, x) 

D 

Here is a 4-step example with A = (1,1, 2, 2). Let 

( 7 ( 77 , x) = Wr [/7(1, a:), h{2, x), h{4:, x), h{5, x),h{n — 2,x)] 
be the corresponding degree n exceptional polynomial. Then, 

(105a; -b 70a;3 -b 84a;3 -b 24:x'^)h{n, x) = 

3 (n - 7)4» - 4). + 5, X) 


16 (77)2(774-3)2 

7 (77-7)2(80- 5777 -b 9772) 


8 (77 - 1)(77 -b 1)2 
105 

h{n -b 3, a;) -I —~ ^)h{n + 1, x) 


4 (n — 1)2 

-b 105(77 — 7 ) 2(77 — 3)2h{n — 1, x) -b 14(77 — 4 ) 3(332 — III 77 -b 9n^)h{n — 3, x) 
+ 84(77 — 7 ) 5 ( 7(77 — 5, x) -b 24(77 — 8 ) 7 ( 7(77 — 7, x) 


7. Closing Remarks 

Connections have already been recognized between the well-established theory 
of bispectrality and the more recent subject of exceptional orthogonal polynomials. 
However, as this note demonstrates, these connections have not been fully utilized. 
Although progress had been made recently on the question of recurrence relations 
for the exceptional orthogonal polynomials, their existence and precise form re¬ 
mained a difficult problem. The use of the bispectral anti-isomorphism reduces the 
problem to very simple (nearly trivial) algebraic manipulation. 
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Certainly, there are more opportunities for a fruitful exchange of ideas between 
these two closely related fields of research that have not yet been realized. In future 
publications, we hope to address the structure of {A/ :/ € 5} as a commutative 
ring of diff’erence operators and its corresponding spectral curve (cf. [12]) as well as 
the implications of the ad-nilpotency of its elements (cf. (3). 
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